The Annals of Applied Probability 

2008, Vol. 18, No. 4, 1326-1350 

DOI: 10.1214/1214/07-AAP487 

© Institute of Mathematical Statistics, 2008 

TWO SUFFICIENT CONDITIONS FOR POISSON 
APPROXIMATIONS IN THE FERROMAGNETIC ISING MODEL 

By David Coupier 

Universite Lille 1 

A d-dimensional ferromagnetic Ising model on a lattice torus is 
considered. As the size of the lattice tends to infinity, two conditions 
ensuring a Poisson approximation for the distribution of the number 
of occurrences in the lattice of any given local configuration are sug- 
gested. The proof builds on the Stein-Chen method. The rate of the 
Poisson approximation and the speed of convergence to it are defined 
and make sense for the model. Thus, the two sufficient conditions 
are traduced in terms of the magnetic field and the pair potential. In 
particular, the Poisson approximation holds even if both potentials 
diverge. 

1. Introduction. Suppose {I a }aer is a finite family of indicator random 
variables, with the properties that the probabilities F(I a = 1) are small and 
that there is not too much dependence between the I a 's. Then, "the law of 
small numbers" says the sum J2 a ^r should have approximately a Pois- 
son distribution. The "birthday problem" and its variants (see Chen [5]), 
the theory of random graphs (see Bollobas [3] for a general reference or the 
famous paper of Erdos and Renyi [11]) and the study of words in long DNA 
sequences (see, e.g., Schbath [18]) are examples in which a law of small num- 
bers takes place. As the situation studied in this paper, these examples can 
be viewed as problems of increasing size (i.e., the cardinality of V tends to 
infinity) in which the sum X^aer la has a Poisson limit. Two methods are 
often used for proving Poisson approximations; the moment method (see 
[3], page 25) and the Stein-Chen method (see Arratia, Goldstein and Gor- 
don [1], Barbour, Hoist and Janson [2] for a very complete reference, or [5] 
for the original paper of Chen). The second one offers two main advantages. 
Only the first two moments need to be computed and a bound of the rate 
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of convergence is obtained. However, the Stein-Chen method requires to re- 
strict our attention to variables which satisfy the FKG inequality [14]. This 
is the case of spins of a ferromagnetic Ising model. 

Let us consider a lattice graph in dimension d > 1, with periodic boundary 
conditions (lattice torus). The vertex set is V n = {0, ... ,n— l} d . The integer 
n will be called the size of the lattice. The edge set, denoted by E n , will be 
specified by defining the set of neighbors V(x) of a given vertex x: 

(1) V(x) = {y^xeV n ,\\y-x\\ p <p}, 

where the substraction is taken componentwise modulo n, || • \\ p stands for 
the L p norm in M. d (1 <p < oo), and p is a fixed parameter. For instance, 
the square lattice is obtained for p = p = 1. Replacing the L\ norm by the 
Lqo norm adds the diagonals. From now on, all operations on vertices will 
be understood modulo n. In particular, each vertex of the lattice has the 
same number of neighbors; we denote by V this number. 

A configuration is a mapping from the vertex set V n to the state space 
{— 1,+1}. Their set is denoted by X n = {— l,+l}" n and called the configu- 
ration set. The Ising model is classically defined as follows (see, e.g., Georgii 
[16] and Malyshev and Minlos [17]). 

Definition 1.1. Let G n = (V n ,E n ) be the undirected graph structure 
with finite vertex set V n and edge set E n . Let a and b be two reals. The 
Ising model with parameters a and b is the probability measure p a b on 
X n = {-1, +1} V " defined by: Vcr € X n , 




where the normalizing constant Z a ^ is such that ^2 a£ x n L i a,b(°~) — 1- 

Following the definition of Malyshev and Minlos ([17], page 2), the mea- 
sure p at b defined above is a Gibbs measure associated to potentials a and b. 
Expectations relative to p a ^ will be denoted by E 0j j,. 

In the classical presentation of statistical physics, the elements of X n are 
spin configurations; each vertex of V n is an atom whose spin is either positive 
or negative. Here, we shall simply talk about positive or negative vertices 
instead of positive or negative spins and we shall merely denote by + and — 
the states +1 and —1. The parameters a and b are, respectively, the magnetic 
field and the pair potential. The model remaining unchanged by swapping 
positive and negative vertices and replacing a by —a, we chose to study only 
negative values of the magnetic field a. Moreover, throughout the paper, the 
pair potential b will be supposed nonnegative. 
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Many laws of small numbers have been already proved for the Ising model. 
They are obtained in two different contexts; the low temperature case (i.e., 
b large enough) and the large magnetic field case. Chazottes and Redig [4] 
have obtained a Poisson law for the number of occurrences of large cylin- 
drical events. Their result concerns the low temperature case and is based 
on an argument of disagreement percolation. Fernandez, Ferrari and Garcia 
[12] have established the asymptotic Poisson distribution of contours in the 
nearest-neighbor Ising model at low temperature and zero magnetic field. 
In the same context, Ferrari and Picco [13] have applied the Stein-Chen 
method to obtain bounds on the total variation distance between the law 
of large contours and a Poisson process. Ganesh et al. [15] have studied the 
Ising model for positive and fixed values of b. Provided the magnetic field a 
tends to infinity, they proved that the distribution of the number of negative 
vertices is approximately Poisson. In Coupier [6], we have generalized this 
result to (deterministic) objects more elaborated than a single vertex, called 
the local configurations; but the parameter b was still fixed. 

Our goal is to extend these Poisson approximations to values of the pair 
potential b not necessarily bounded. In particular, we establish a law of small 
numbers in the case where both magnetic field a and pair potential b diverge. 
Moreover, compared with the previous articles, the results of this paper offer 
two advantages. First, the roles of the rate of the Poisson approximation 
and the speed of convergence to it are conceptually clearer. Second, the 
method of the proof is more elementary and simple than those of the previous 
references; it does not require the need of auxiliary birth-and-death processes 
as in [12] and [13], or of activity expansions as in [15]. 

Pattern recognition and image denoising problems motivate the study 
realized in this paper. Indeed, Desolneux, Moisan and Morel (see [9] and 
[10]) describe and detect the geometric structures of an image through the 
notion of "meaningful events." Their method is based on the link between 
the perception threshold of a given visual structure and its probability of 
appearance in a random image. In Coupier, Desolneux and Ycart [7], a 
Poisson approximation for the probability of appearance of any local pattern 
is computed in the case of independent pixels. This result leads to a denoising 
algorithm for gray-level images. However, the structures created by the noise 
are generally more irregular (with a large perimeter) than the natural ones 
of an image. This remark justifies our interest for Poisson approximations 
in probabilistic models with dependent pixels, as the Ising model. 

We are interested in the occurrences in the graph G n of a fixed local 
configuration i] (see Section 2 for a precise definition and Figure 1 for an 
example). Such a configuration is called "local" in the sense that the vertex 
set on which it is defined is fixed and does not depend on n. Its number of 
occurrences in G n is denoted by X n (rj). A local configuration rj is determined 
by its set of positive vertices V+iri) whose cardinality and perimeter are 
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respectively denoted by k{rj) and 7(7/). As the size n of the lattice tends to 
infinity, the potentials a = a(n) and b = b{n) depend on n. The case where 
a(n) tends to —00 corresponds to rare positive vertices among a majority 
of negative ones. In order to simplify formulas, the Gibbs measure /J a (n),&(n) 
will be merely denoted by fj, a ,b- 

A natural idea consists in regarding both potentials a(n) and b(n) through 
the same quantity: the weight of the local configuration 77 

W n ( V ) =exp(2a(n)A;(7?) - 26(71)7(77)). 

Actually, the weight W n {rj) represents the probabilistic cost associated to a 
given occurrence of 77 and, consequently, the product n d W n (rj) (where n d is 
the cardinality of V n ) corresponds to the expected number of occurrences of 
77 in the whole graph G n . Hence, in order to obtain a Poisson approximation 
for the random variable X n (rj) it is needed that the weight W n (rj) tends to 
at the rate n~ d . Therefore, throughout this paper, the potentials a[n) and 
b{n) will satisfy the homogeneity relation 

(3) n d W n {r,) = A, 

where A is a positive constant. Two other parameters naturally take place in 
our study: the probability gap A n (ry) and the maximality probability Q n (r]). 
When the first one tends to 0, the vertices surrounding a given occurrence 
of 77 are all negative. A null limit for the second one means if somewhere 
in the graph the vertices corresponding to the set V+ (77) are positive, then 
locally they are the only ones. The fact that the probability gap A n (r/) 
and the maximality probability 6^(77) tend to is denoted by (HI). We 
denote by (H2) the following (interpreted) hypothesis: "the probability for a 
given configuration £ G { — ,+} v , V C V n , to occur is bounded by its weight 
W n {Q. n Our main result (Theorem 3.1) describes the asymptotic behavior 
of the number of occurrences of 77 in the lattice. Added to the necessary 
condition (3), the two hypotheses (HI) and (H2) imply the total variation 
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Fig. 1. A local configuration r\ with k(n) = |V+(r;)| = 10 positive vertices and a perimeter 
7(77) equal to 58, in dimension d = 2 and on a ball of radius r — 2 (with p = 1 and relative 
to the Lao norm). 
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distance between the distribution of the random variable X n {rj) and the 
Poisson distribution with parameter A satisfies 

d T y(C(X n (7 ] )),V(X)) = 0(m a x{A n (r ] ),Q n ( V )}). 

The proof is based on the Stein-Chen method. Lemma 4.1 reduces the proof 
to a sum of increasing random indicators X n {rj). Since the Gibbs measure 
/j,a t b defined in (2) satisfies the FKG inequality (because the pair potential 
b is positive), this method is applied to the random variable X n (rj) and 
produces Lemma 4.5. Then, bounds on the first two moments of X n (rj) 
(Lemmas 4.2 and 4.6) allow to conclude. 

Finally, Theorem 3.1 is completed by a second result (Proposition 3.2) 
which traduces the two hypotheses (HI) and (H2) in terms of magnetic field 
a{n) and pair potential b(n). 

The paper is organized as follows. The notion of local configuration 77 is 
defined in Section 2. Its number of positive vertices k(rj), its perimeter 7(7/) 
and its weight W n {rj) are also introduced. The use of the local energy of 7/ 
(Definition 2.1) allows to write the conditional probability for rj to occur in 
the graph in terms of weights (Lemma 2.2). Lemma 2.4 is the main result 
of Section 2. It provides lower and upper bounds for this conditional proba- 
bility on which the estimates of the first two moments of A^ n (?7) are based. 
In Section 3, Theorem 3.1 and hypotheses (HI) and (H2) are introduced 
and discussed. Thus, Proposition 3.2 is illustrated by two examples. Finally, 
Section 4, in which the sum X n {rf) is defined, is devoted to the proof of 
Theorem 3.1. 

2. Conditional probability of a local configuration. Let us start with 
some notation and definitions. Given £ G X n = { — ,+} Vn and V C V n , we 
denote by (v the natural projection of £ over {— ,+} v ■ If U and V are two 
disjoint subsets of V n , then QuC'v ls the configuration on U U V which is 
equal to £ on U and £' on V. Let us denote by 5V the neighborhood of V 
[corresponding to (1)]: 

SV = {y G V n \ V 3x G V, {x,y} G E n }, 

and by V the union of the two disjoint sets V and SV. Moreover, | V | denotes 
the cardinality of V and J~(V) the cr-algebra generated by the configurations 
of{-,+} V . ' 

As usual, the graph distance dist is defined as the minimal length of a 
path between two vertices. We shall denote by B(x,r) the ball of center x 
and radius r: 

B(x, r) = {y£ V n ; dist (a;, y) < r}. 

In the case of balls, B(x,r) = B(x,r + 1). In order to avoid unpleasant sit- 
uations, like self-overlapping balls, we will always assume that n > 2pr. If n 
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and n' are both larger than 2pr, the balls B(x,r) in G n and G n > are isomor- 
phic. Two properties of the balls B(x,r) will be crucial in what follows. The 
first one is that two balls with the same radius are translates of each other: 

B(x + y,r) =y + B(x,r). 

The second one is that for n > 2pr, the cardinality of B(x,r) depends only 
on r and neither on x nor on n: it will be denoted by (3{r). The same is 
true for the number of edges {y,z} E E n with y,z £ B(x,r), which will be 
denoted by a(r). 

Let r be a positive integer, and consider a fixed ball with radius r, say 
B(0, r). We denote by C r = {—, -\-} B (°> r ) the set of configurations on that ball. 
Elements of C r will be called local configurations of radius r, or merely local 
configurations whenever the radius r will be fixed. Of course, there exists 
only a finite number of such configurations (precisely 2^( r )). See Figure 1 
for an example. Throughout this paper, the radius r will be constant, that 
is, it will not depend on the size n. Hence, defining local configurations on 
balls of radius r serves only to ensure that studied objects are "local." In 
what follows, rj, r( will denote local configurations of radius r. 

A local configuration 77 6 C r is determined by its subset V + (rj) C B(0,r) 
of positive vertices: 

V + (ri)={x£B(0,r), V (x)=+}. 

The cardinality of this set will be denoted by k(rj) and its complementary 
set in B(0,r), that is, the set of negative vertices of 77, by V-(rj). Moreover, 
the geometry (in the sense of the graph structure) of the set V+{rj) needs to 
be described. Let us define the perimeter 7(77) of the local configuration rj 
by the formula 

7(7?) = V\V+( V )\ - 2\{{x,y} E V + ( V ) x V + ( V ), {x,y} E E n }\, 

where V is the number of neighbors of a vertex. In other words, 7(7/) counts 
the pairs of neighboring vertices x and y of 5(0, r) having opposite spins 
(under 77) and those such that x E -8(0, r), y E 6B(0, r) and rj{x) = +. Finally, 
we denote by W n {rj) and call the weight of the local configuration r\ the 
following quantity: 

W n (rj) = exp(2a(n)k(r]) — 26(71)7(77)). 

Since a(n) < and b(n) > 0, the weight W n (r]) satisfies < W n (-q) < 1. That 
of the local configuration having only negative vertices, denoted by 77 and 
called the null local configuration, is equal to 1. If 77 7^ 77 , then A; (77) > 1 and 
7(77) > V. It follows that 

W n (77) < exp(2a(n) - 26(n)V). 
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Actually, the weight W n (rj) represents the probabilistic cost associated to 
the presence of r\ on a given ball. Lemma 2.4 will give a rigorous sense to 
this idea. 

Remark that the notation &(■), and W n (-) can be naturally extended 
to any configuration £ € { — , +} , V C V n . Furthermore, the configuration 
of {— , +} v having no positive vertex will be said to be null and denoted by 

C°- 

Let j] 6 C r . For each vertex x E V n , denote by n x the translation of r\ onto 
the ball B(x,r) (up to periodic boundary conditions): 

Vy€V n dist(0,y)<r rj x (x + y) = rj(y). 

In particular, V+(i] x ) = x + V+(r]). So, r\ and r] x have the same number 
of positive vertices and the same perimeter. So do their weights. Let us 
denote by the indicator function defined on X n as follows: I^(cr) is 1 if 
the restriction of the configuration a £ X n to the ball B(x,r) is rj x and 
otherwise. Finally, let us define the random variable X n {rf) which counts the 
number of copies of the local configuration r/ in G n : 

Xn(v) = E ^ 

xGV„ 

Due to periodicity, this sum consists of n d indicator functions IJJ, which 
have the same distribution. Hence, in order to understand the behavior of 
the random variable X n {rf) it suffices to study that of one of the random 
indicators I^- This is the reason why the rest of this section is devoted to 
the study of the probability for r/ (or rj x ) to occur on a ball B(x,r) knowing 
what happens on its neighborhood: 

(4) Ha,b{I x i = l WsB{x,r))- 

Besides, let us remark that the expected number of occurrences of the local 
configuration r\ in the graph can be expressed as the following expectation: 

K ajb [X n ( V ))=E atb [n d I2] 

= E a}b [n d v a , b (l2 = l\T(6B))}, 

where x is a given vertex and B denotes merely the ball B(x,r). Here, 
^i a ^{I^. = 1\J-(5B)) represents a J 7 (5B)-measurable random variable and, 
for any configuration a of {— ,+} SB , /J, ajb (I^ = l\J r (5B))(a) is equal to the 
conditional probability fj, ajb (I^ = l\a). 

Lemma 2.4 provides inequalities between the conditional probability (4) 
and the weight W n (jj). A way to link these two quantities consists in using 
the local energy of rj. Let us start with the following definition. 
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Definition 2.1. Let V cV n and C G The /oca/ eneroy i/ y (C) 

of the configuration £ on the set V is defined by 

H v (C)=a(n)Y J C(y)+b(n) £ C(l/)C(*), 
yeV {y,z}€E n 
(yeV)v(zeV) 

where (y £ V) V (z S V) means at least one of the two vertices y and z 
belongs to V (the other might belong to its neighborhood SV). 

Let us fix a vertex x and denote by B the ball B(x,r). On the one hand, 
remark that for any a E { — ,+} SB the local energy H B (r] x cr) on B of the 
configuration which is equal to r] x on B and er on SB gives an explicit 
formula of the conditional probability /U a &(Ij[? = 1|<t): 

/-x , /ti, 1 1 \ exp(ij' i? (? ?x o-)) 

I 6 ) ^,6(4 = 1 F) = ^ ,TT B ( , w 

On the other hand, the exponential of the local energy of rj can be expressed 
as a function of the weight W n (rj). Indeed, H B (r] x a) is equal to 

(7) a(n)(2k(ri)-P(r)) + b(n)( £ Vx (y)ri x (z) + ]T %(y)a(z)J. 

\{y,2}e£; n {j/,^}ge„ / 

y,z£B y£B,zeSB 

Let {y,z} be an edge such that y £ B and ^(y) = +. Assume z £ B. If 
TfaOz) — — j the edge {y,z} is counted by the perimeter 7(77), whereas if 
Vx( z ) = + , it is not. In the case where z £ B (i.e., z G SB), the edge {y, z} is 
consistently counted by 7(7?). In other words, the following equality holds: 

(8) 27(77) =a(r)- J2 Vx(y)rix(z) + £ ( 1 + Vx(y))- 

{y,z}GE n {y,z}&E n 
y,z£B y£B,ze5B 

So, the perimeter 7(7/) can be inserted in the expression of H B {rj x a) given 
by (7) and we get 

exp(H B ( Vx a)) 

(9) =W n (r7)exp(-a(n)/3(r) 



+ 6(n)(a(r)+ £ l + r) x (y)(l + a(z))\). 

V {y,z}£E n ) ) 

y£B,z£8B 

Remark that the previous identity remains valid replacing the ball B with 
any given set V and the configuration r\ x a with any given £ € { — , +} • 
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Using the previous relations, an alternate expression of the conditional 
probability iJi a fi{I x = l|a) is obtained. 

Lemma 2.2. Let x G V n and B = B(x,r). Then, for any rj G C r and 

W n (rj x a SB ) 



(10) fi a , b (I^ = l\a SB ) 



T, v '£C r W n(r]' x O-5B)' 



Proof. Let rj be a local configuration of radius r and a G {— 
Let us denote by D the set 5B U 55: a vertex y belonging to D satisfies 
r < dist(x, y) < r + 2. Adding terms which only depend on the configuration 
a to the right-hand side of (6) does not change the equality: 

„ /r ,_ lU ) exp(H B (r] x a D )) 

(11) Va,b{lx = M a SB) ~ 



E v '£C r ex P( HB (Vx a D)) 

Let 7/ G C r . Applying relation (9) to the configuration rj' x ar>, we get 
exp(H*( V > x a D )) 

= W n (rj x a SB ) exp ^-a(n)(3(r + 1) 

+ b(n)(a(r + l)+ £ 1 + a D (y)(l + a D {z)) ) ) . 
y£B,ze6B 

Finally, it suffices to remark that the latter exponential does not depend on 
rj' x and to simplify the right-hand side of (11). □ 

Let V and V be two disjoint subsets of vertices. The following relations: 

fc(CO = fc(C) + fc(0 and 7(CO< 7(0 +7(0 

are true whatever the configurations £ G { — ,+1^ and £' G { — ,+} V ■ The 
connection between ( G {— ,+} v and £' G {— , denoted by conn(£,£'), 
is defined by 

conn(C,C0 = \{{y,z} €E n ,yeV,zeV and ((y) = ('(z) = +}|. 

This quantity allows to link the perimeters of the configurations C an d 
(' together: 

7(CO + 2conn(C,0 = 7(C)+7(0, 
and as a consequence their weights: 

(12) W n (((') exp(-46(?i) conn(C, O) = W n (()W n (('). 
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In particular, if the connection conn(£,£') is null, then the weight W n (CC') 
is equal to the product W n {QW n {C,') . This is the case when V CiV = . 

Let B = B(x,r) and rj be a local configuration of radius r. Assume the 
set of positive vertices of rj satisfies V+(rj) C 5(0, r — 1); such a local con- 
figuration is said to be clean. Then, the connection coim(r] x , cr) is null and 
W n (rj x a) is equal to the product W n (rj)W n (a). This is the case also for the 
null configuration cr°; W n (r/ X a ) = W n (r])W n (a ) = W n {rf). 

Let us define the probability gap of w and denote by A n (7/) the quantity 
below: 

a / \ W n (r] x a) 
A„ m = max . . , 

Since the weights of a configuration and its translates are equal, the prob- 
ability gap A n (77) does not depend on the ball B = B(x,r) (nor on 5B). It 
only depends on the local configuration r\ and on the potentials a(n) and 
b(n). Moreover, A n ,(r/) satisfies the following inequalities: 

Lemma 2.3. Let r]^C r . Then, the probability gap A n (?y) satisfies 
(13) exp(2a(n) - 2b(n)V) < A n (r?) < exp(2a(n)). 

In particular, for all r( sC r , rj ^ rp : 

W n (rf) < A n (r]). 

Proof. Let rj £ C r and as usual denote by B the ball B(x,r). Since 
o(n) < and b(n) > 0, the weight W n (r] x a) is a decreasing function of pa- 
rameters k(r] x a) and j(r/ x a). If a and a' are two configurations of { — , +} SB , 
different from cr°, such that V + (a) C V + {a') and k(a) = k(a') — 1, then the 
perimeter 7(?7 x cr) is not necessarily smaller than 7(77 x cr'). However, the fol- 
lowing statement based on the convexity of balls on which the local configu- 
rations are defined is true: for all configuration a G {— , +} SB , a 7^ a , there 
exists a configuration cr m i n (a) € { — ,+} SB satisfying 

k(o- min (a)) = 1, V + (a min (a)) C V+(cr) and 7(7? a; (T in i n (cT)) < 7(^0-). 

As a consequence, the maximum of W n (r] x o~) is obtained among the configu- 
rations a € {— , +} SB having only one positive vertex; k(a) = 1 and 7(a) = V. 
For such a configuration cr, 

W n (rj x a) 



W n (rj) 



= exp(2a(n) - 2b(n)(j(rj x a) - 7(7?))) 

= exp(2a(n) — 2b(n)(V — 2 conn (r) x , cr))) 
>exp(2a(n) -26(n)V). 
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The lower bound of (13) follows. The upper bound is also a consequence of 
the convexity of the ball B: the perimeter 7(77^0"), a 7^ a , is necessarily as 
large as j(r) x a°). In other words, the difference 7(77^0") — 7(7/) is positive. 

Finally, we have already seen that the weight of a configuration ?/ £ C r , 
different from r/°, is bounded by exp(2a(n) — 2b(n)V). So, it is bounded by 
A n ( v ). □ 

For example, the probability gap of a clean local configuration is equal 
to the lower bound given by (13), that is, exp(2a(n) — 2b(n)V). Indeed, if r\ 
is clean, the connection conn(r/ x , •) is null and the equalities of the previous 
proof allow us to conclude. 

Lemma 2.4 is the main result of this section. It provides lower and up- 
per bounds for the conditional probability n a j,(LJj. = l|cr) depending on the 
weight W n {rj) and the probability gap A n (ry): 

Lemma 2.4. Let r\ be a local configuration of radius r and let x be a 
vertex. Let us denote by B the ball B(x,r). Then, for all configuration a £ 

{-,+}"*, 



(14) /i a , 6 (I£ = \\a) > W n ( V )t a=a o(l - \C r \A n (n)) 
and 

(15) ^ b (P x =l\a)< W n ( V ) (l + t a ^ . 

Under the hypotheses of Theorem 3.1, the probability gap A n (r?) tends to 
and the probability of the indicator l cr=(T o tends to 1. Henceforth, relations 
(14) and (15) mean the weight W n {rf) constitutes a good estimate for the 
probability /i a , 6 (I^ = 1). 

Before proving Lemma 2.4, let us explain why the results of this section 
do not hold when the pair potential is negative. If b(n) < 0, the condition 
a{n) — b{n)V < (which ensures that the weight of a local configuration is 
smaller than 1) implies the weight W n (rj) and the probability gap A n ,(ry) are 
both bounded by exp(2a(n) — 2b{n)V). In particular, we lose the inequality 
W n {r)') < A n (rj) of Lemma 2.3 and therefore the lower bound of Lemma 2.4. 

Proof OF Lemma 2.4. Let 77 £ C r and x £ V n . Let a be a configuration 
of { — , +} SB . If a = a°, then Lemma 2.2 implies that 



/v,(i2 = iK 



W n (rj) 
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W n (rj) 



<w n ( v ), 

whereas, if a ^ a , 

W n (r] x a) 1 



lia, b (l2 = l\a) = W n (ri)- 



< W n (r ? )A n (r ? ) 

< W n ( V ) 



W n (v) E v >eC r W n (r,> x a) 
1 



W^) 

An.(??) 



W n ((T)' 

The upper bound (15) is deduced from the above inequalities. The lower 
bound is also based on Lemma 2.2. Indeed, 

(tv 1 1 \ ^ -n W n (r] x a°) 



WM 



Now, the weight of a local configuration 77' 7^ 77° satisfies W n (?/) < A n (r/) 
(see Lemma 2.3). Hence, we can write 

/i a ,fe(^ = 1|<T) >l CT=(T o- 



l + |C r |A„(7?)' 

Finally, the lower bound (14) is obtained using the classical inequality 

Vu>-1 — - — >l-u. n 
1 +u ~ u 

3. Poisson approximation. This section is devoted to the main result of 
the paper, that is, Theorem 3.1. Let us introduce two hypotheses (HI) and 
(H2) that will be sufficient to prove Theorem 3.1. 

For any given local configuration rj, let us define the subset C r {rf) of C r by 

C r ( V ) = W£C r ,V + (i 1 >)DV + ( V )}. 

Each element of C*{rf) = C r (jj) \ {r/} has at least k(rj) + 1 positive vertices. 
Thus, we denote by O n (ry) and we call the maximality probability the fol- 
lowing quantity: 

e („\- max W »W) 
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The necessary condition n d W n (r]) = A, for some constant A > 0, forces the 
weight W n {rf) to tend to and therefore the local configuration rj to be 
different from rp. Thus, using the inequality 7(77) < Vk(r]), we deduce that 
the probabilistic cost associated to a single positive vertex, that is, 

exp(2a(n) -2b(n)V), 

is bounded by W n {ri) l / k ^ and tends to as n tends to infinity. Actually, 
Theorem 3.1 requires stronger hypotheses. First, if the vertices correspond- 
ing to x + V+(j]) are positive, then the other ones belonging to B(x,r) and 
those belonging to the neighborhood SB(x,r) must be negative. This fact 
results in 

(HI): lim max{A n (r ? ),9 n (7 ? )} = 0. 

n— >+oo 

For V C V n and £ G { — , +} v , the indicator function Iy is defined as follows: 
Iy{<j) is 1 if the restriction of the configuration a G X n to V is C an d 
otherwise. Theorem 3.1 needs a control of the probability of the event Iy = l. 
Let (H2) be the following hypothesis: 

W C V n , 3N(V) G N, 3C = C(\V\) > 0, VC G {-, +} v , 

{ j ' Vn>N(V),n a , b (I v = l)<CW n (0- 

Here, C = C(\V\) means that C depends on the set V only through its 
cardinality. 

Before introducing our main result, let us recall some classical notations. 
If /i and v are two probability distributions, the total variation distance 
between \i and v is 

cZ T v (m, v) = sup I n(A) - v(A) I, 
A 

where the supremum is taken over all measurable sets. The probability dis- 
tribution of a random variable X is denoted by C{X) and that of Poisson 
distribution with parameter A by "P(A). Moreover, if f(n) and g{n) are two 
positive functions, notation f(n) = 0(g{n)) means that there exists a con- 
stant C > such that, for all n, f(n) < Cg(n). 

Theorem 3.1 states that the asymptotic behavior of X n {rf) is Poissonian 
and defines the speed of convergence in terms of total variation distance. 

Theorem 3.1. Let r/ be a local configuration of radius r. Assume that 
the magnetic field a{n) is negative, the pair potential b{n) is nonnegative 
and they satisfy 



n d W n ( V ) = A, 
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for some constant A > 0. Furthermore, if the hypotheses (HI) and (H2) are 
satisfied, then the total variation distance between the distribution of X n {rf) 
and the Poisson distribution with parameter X is such that: 

ckv(C(X n (r ) )),V(\)) = 0(m a x{A n (r ] ),e n (r))}). 

Now, some comments are needed. First, this result generalizes the Pois- 
son approximation given in [6] (Theorem 1.3). When the pair potential 
b(n) = b > is fixed, the product n d exp(2a(n)/c(r/)) becomes constant. De- 
note by A' this quantity. Then, the limit distribution for X n {rf) is the Pois- 
son distribution with parameter A = A' exp(— 2b{n)^(rj)). Furthermore, the 
probability gap A n (ry) and the maximality probability Q n (r]) are of order 
exp(2a(n)) and consequently, 

d TY (£(X n ( V )),V(X'eM-nnh(v)))) = 0{n~ d ' k ^). 

We believe that max{A n (r/), Q n (j])} is the real speed at which the total 
variation distance between C{X n {rj)) and V{X) tends to zero. Indeed, it 
seems to be true for the upper bound given by Lemma 4.5 (for more details, 
see Chapter 3 of [2]). Moreover, this has been proved by Ganesh et al. [15] in 
the case where the local configuration rj represents a single positive vertex 
[with k{rj) = 1, 7(7?) = V] and the pair potential is fixed. 

Under the hypotheses (HI) and (H2), the expectation of X n (jj) satisfies 

n d W n (r,)(l - 0(A n ( V ))) < E a , b [X n (r})] < n d W n (n)(l + 0(A n ( V ))) 

[see relations (20) and (21) of the next section]. In particular, if the product 
n d W n (ij) tends to (resp. +00), the same is true for E a ^[X n (rj)]. Actually, 
a better result can be easily deduced from Theorem 3.1: if the product 
n d W n (rj) tends to (resp. +00), then the probability ii a ^{X n {n) > 0) tends 
to (resp. 1). Roughly speaking, if W n {ifj is small compared to n~ d , then 
asymptotically, there is no occurrence of rj in G n . If W n (jj) is large compared 
to n~ d ', then at least one occurrence of r\ can be found in the graph, with 
probability tending to 1. Using the vocabulary of the random graph theory, 
this means that the quantity n~ d is the threshold function for W n (r]) of the 
property X n {rj) > 0. A straight proof of this statement can be found in [8]. A 
result similar to (actually weaker than) Theorem 3.1 can be proved without 
using the Stein-Chen method. Under the same hypotheses, the distribution 
of X n {rj) converges weakly to V(X), as n tends to infinity. The proof of this 
result is based on the moment method (see [3], page 25) and Lemma 2.4; 
see the ideas developed in the proof of Lemma 4.6 for details. 

Replacing Proposition 4.4 with the Stein-Chen formulation of [1], it seems 
the estimates given by Lemma 2.4 should apply to systems formed by "ani- 
mals" interacting only by volume exclusions (as in [12] and [13]). 
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Finally, as has been suggested by an anonymous referee, Theorem 3.1 
can be extended to some nonlocal configurations whose cardinalities diverge 
suitably with the size of the lattice. 

The rest of this section is devoted to translating the hypotheses of Theo- 
rem 3.1 in terms of magnetic field a(n) and pair potential b(n). 

Proposition 3.2. Let n be a local configuration of radius r. Assume 
that the magnetic field a(n) is negative, the pair potential b(n) is nonnegative 
and they satisfy 

n d W n (? ] )=X, 

for some constant A > 0. Then the two conditions a(n) + 2Vb(n) < 0, for n 
large enough, and 

lim a(n) + Vb(n) = — oo 

n— >+oo 

imply respectively the two hypotheses (H2) and (HI). 

The limit a(n) + Vb{n) — > — oo forces the magnetic field a(n) to tend to 
— oo. Proposition 3.2 and Theorem 3.1 state that a(n) — > — oo is the main 
condition bearing on the magnetic field a(n) in order to expect a Poisson 
approximation. Rare positive vertices among a majority of negative ones 
constitutes the global context in which our Poisson approximations take 
place. Conditions bearing on the pair potential b(n) are much larger. Here are 
two examples of couples (a(n),b(n)) satisfying the conditions of Proposition 
3.2 and so the hypotheses of Theorem 3.1. 

Example 3.3. Assume the potentials a(n) < and b(n) > are such 
that n d W n (i]) = A, for some A > 0, and 

bin) 

lim a{n) = — oo and lim — — — - = 0. 

n-»+oo n-»+oo a(n) 

Then, conditions of Proposition 3.2 are satisfied. Indeed, the quantity 

a(n) 



a(n) + Vb(n) = a{n)[ 1 + V 



tends to — oo, so does a(n) + 2Vb{n). The pair potential b{n) is allowed to 
tend to +oo [but slower than |a(n)|] or to be bounded. 

Example 3.4. Assume the potentials a(n) and b(n) satisfy 
a(n)= 2Mr ? )+7(r/)/V l0g (^' 



h{n ) = = - \oJ- 
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where A is a positive constant. The rate at which the pair potential b(n) 
tends to infinity is the same as |a(n)|. It is easy to check that the magnetic 
field a(n) is negative, b(n) is nonnegative, the product n ,W n (rj) is equal to 
A, the quantity a{n) + 2Vb(n) is null and a(n) + Vb(n) = a(n)/2 tends to 
— oo. In other words, conditions of Proposition 3.2 are satisfied. 

Let us end this section by proving Proposition 3.2. 

Proof of Proposition 3.2. Let r\ be a local configuration of radius r. 
Assume that potentials a(n) and b(n) satisfy conditions of Proposition 3.2. 

First, let us prove (HI) is satisfied. The study of the probability gap has 
been already done; A n (ry) < exp(2a(n)). Hence, the limit a(n) + Vb(n) — > 
— oo forces the magnetic field a(n) to tend to — oo and A n (r/) to 0. Now, let 
us bound the maximality probability & n (rj). Let rf G C*(rj). Although the 
set of positive vertices of rj' contains that of rj, the perimeter 7(77') is not 
necessarily larger than 7(77): roughly, V+(r?) may have holes. However, the 
inequality 

7(r?')>7(r?)-(^0-^))V 
holds. Hence, the ratio W n (rf) divided by W n (rj) is bounded: 

-^U- = exp(2a(n)(fc(r? / ) - fcfa)) - 2b(n)( 1 (r j ') - 7 (r?))) 

< exp((2a(n) + 2Vb{n)){k{r]') - k{r)))). 
Then the limit a(n) + Vb(n) — > —00 implies the maximality gap satisfies 

e„<„ )= , mx ^W) 

v'ecf( v ) W n (rj) 

< exp(2a(n) + 2Vb(n)) 

and tends to as n tends to infinity. 

It remains to prove that (H2) holds. Let V be a set of vertices and let £ 
be a configuration of {— , +} ■ Let us start with writing 

/v fe (4 = i)= £ /x a , 6 (4=j,v=i). 

Let a and to be two configurations respectively defined on 5V and SV. 
Thanks to Lemma 2.2 and relation (12), we get 



< W n (CaLo) 



W n {0j) 
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< W n ((a) exp(46(n) conn(C<7, u)) 

< W n (C<7) exp(46(n) conn(cr, u)). 

Moreover, the weight of the configuration Qa can be expressed as 

W n ((a) = W n (0 exp(2a(n)k(a) - 26(n)( 7 (C C r) - 7 (C))) 

= W n (C)exp(2a(n)A;(cr) - 26(n)(7(cx) - 2conn(C,cr))). 

Combining the previous relations, we deduce the following inequality: 

/v 6 (4 CT =i|/^=i) 

< W n (()exp(2a(n)k(a) - 2b(n)(j(a) - 2conn(C,cr) - 2conn(cr, w))). 

A way to obtain an upper bound for the above conditional probability which 
does not depend on u consists in using the inequality 

conn(£,0") + conn(<r, uj) < Vk(a). 

As a consequence, 

^ b (l£ = l\Tfo = 1) < W n (C) exp((2a(n) + AVb(n))k(a) - 2b(n)y(*)) 
<W n (Q, 

whenever a(n) + 2Vb{n) is negative. Then, the probability ^L a ^(I^ = 1) is 
bounded by the weight W n (C) an d finally we get 

/ia, b (4 = l)<2 |5V Vn(C)- 

Observe that the constant C(V) = 2 V '^I > 2^ is suitable and depends on 
the set V only through its cardinality. □ 

4. Proof of Theorem 3.1. This section is devoted to the proof of The- 
orem 3.1, whose layout is essentially the same as the one of the proof of 
Theorem 1.3 of [6] [i.e., in the case b(n) = b > 0]. Let ij be a local configu- 
ration of radius r. Throughout this section, we shall suppose the magnetic 
field a(n) is negative, the pair potential b{n) is nonnegative and they satisfy 

n d W n {ri) = A, 

for some constant A > 0. 

For all vertex x, the indicator random variable 7^ is defined by 

T= V rf 

± X / j X 1 

r)'eC r (ri) 
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where C r (r/) is formed by the local configurations of radius r whose set of 
positive vertices contains that of rj. Thus, we denote by X n (rj) the sum of 
these indicators over x S V n : 

X n (ri) = £ T x 

xev„ 

(16) 

= X n ( V )+ ]T X n (rf), 
v'ec*( v ) 

whose expectation E aj fe[X n (ry)] will be simply denoted by A n . Finally, in 
order to simplify formulas, the maximum of the probability gap A n (?y) and 
the maximality probability Qniv) wm be merely denoted by M n (rj): 

M n {rj) = max{A n (r/), @ n (v)}- 

The proof of Theorem 3.1 is organized as follows. The total variation distance 
between £(X n (r])) and 'P(A) is bounded by 

d T v(>C(X rl (77)) )J C(X n (7 ? )))+dTv(>C(X n (7 7 )),P(A n )) + dTv(^(An),P(A)). 

We are going to prove that each term of the above sum is of order 0(M n (n)). 
The first and the last ones are respectively dealt with using Lemmas 4.1 and 
4.3. Applied to the family of indicators {1%, x G V n }, the Stein-Chen method 
gives an upper bound for the second term (Lemma 4.5). Finally, Lemma 4.6 
implies that this upper bound is a C(M n (r/)). 

As the maximality probability @ n (i]) tends to 0, the occurrences of local 
configurations of C*(r]) = C r {rf) \ {r/} have vanishing probability. Hence, the 
random variables X n (rj) and X n {j]) (resp. their expectations E[X n (r/)] and 
X n ) will be asymptotically equal. 

Lemma 4.1. The total variation distance between the distributions of 
X n (n) and X Tl {rj) satisfies 

(17) d TY (£(X n ( v )),£(X n (v)))=0(M n ( v )). 
Furthermore, the following inequalities hold: 

(18) K,b[ X n(v)] <\n< KA X n(v)} + K (v) I AC( J B)9 n ( V ), 

where C(B) denotes the constant of the hypothesis (H2) corresponding to a 
ball of radius r. 

Proof. The total variation distance between two probability distribu- 
tions \i and v can be written as 

drv(fi, u) = inf{P(X ^ Y),£(X) = fi and £(Y) = u}. 
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Using this characterization and the identity (16), it follows that 
ckv(C{X n (n)),£(X n (r)))) < n a , b {X n {ri)^X n {r,)) 

< Ha,b(X n (v) > X n {r,)) 



\ V 'eCf(r,) J 



The above sum is an integer-valued variable. So, its probability of being 
positive is bounded by its expectation. Hence, 

<hv(C(X n (n)),C(X n (ri)))< E KA x n(r)')}. 

rfeC*{r,) 

Let rf sC*(ry). The hypothesis (H2) allows to control the expectation of rj' 
with the maximafity probability Q n (rj): 

E a>6 [X n (V)]=nV6(^' = l) 
< n d C(B)W n (rf) 

W n (r{) 
<XC(B)e n ( v ). 

Finally, we get 

d TY (C(X n ( v )),C(X n ( v ))) < \e*(r,)\XC(B)@ n (rj) 

and (17) follows. Relation (18) is an immediate consequence of the previous 
inequalities and the identity (16). □ 

Lemmas 2.4 and 4.1 state that the expectation of X n (r]), that is, A n , is 
close to that of X n (rj), which is itself close to A. 

Lemma 4.2. There exists a constant K(r) > such that 
(19) A(l - K(r)M n (n)) < A„ < A(l + K(r)M n ( V )). 

The constant K(r) depends on the radius r and on the parameters of the 
model (i.e., p, p and the dimension d). 

PROOF of Lemma 4.2. Let xeV n and denote by B the ball B(x,r). 
First, note that the indicator t^^o occurring in Lemma 2.4 can be expressed 
as the following sum: 

E ^B' 
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where cr° represents the null configuration of {— , +} SB . So, the upper bound 
(15) implies 

<\(l + A n ( v ) £ E ^]\ 



Let C{5B) be the constant of the hypothesis (H2) corresponding to the 
neighborhood of a ball of radius r. Then, 

KA^b) < C(5B)W n (a). 

Hence, we obtain an upper bound for E 0j {,[X n (7/)]: 

(20) VaAMv)] < A(l + A n ( v )2 SB C(5B)). 

Let us deal with the lower bound. It has been already said (Lemma 2.3) 
that the weight of a configuration a different from a is smaller than the 
probability gap A n (rj). As a consequence, 



E 



a.b 



< £ C(£B)W n (<x) 
<2 SB C{SB)A n (rj). 



Next, we deduce from the lower bound (14): 

E a:b [X n ( V )} =E , 6 [n ( Va,&(^ = 11^(55))] 

>A(l-|C r |A n (? ? ))^l-E a 

(21) 



2 T SB 



> A(l - |C r |A n ( ?? ))(l - 2 SB C(6B)A n (rj)) 
>X(l-A n ( v )(\C r \ + 2 SB C(5B))). 

Finally, combining the inequalities between the expectation of X n [rj) and 
A n given by relation (18) and the inequalities between the expectation of 
X n {rj) and A given by (20) and (21), we conclude with 

K(r) = max{|C r | + 2 SB C(5B), \e*(rj)\C{B) + 2 m C{5B)}. □ 

The total variation distance between two probability distributions on the 
set of integers can be expressed as 

^tv(M) v ) = \ Yl \v( m ) ~ v{m)\. 

m>l 



POISSON APPROXIMATIONS FOR THE ISING MODEL 



21 



Using this characterization, relation (19) between A n and A and hypothesis 
(HI), it is possible to bound the total variation distance between 7'(A n ) and 
■p(A). This result is very similar to Lemma 4.2 of [6]. So, it will not be proved 
here. 

Lemma 4.3. The total variation distance between the Poisson distribu- 
tions with parameters A n and A satisfies 

d TY (V(X n ),V(X)) = 0(M n ( v )). 

There remains to bound the total variation distance between C(X n (rj)) 
and V(X n ). This is based on the Stein-Chen method and particularly on 
Corollary 2.C.4, page 26 of [2] which is described below (Proposition 4.4). 
Let {Ij}j g / be a family of random indicators with expectations 7Tj. Let us 
denote 

Z = Ii and 9 = 7Tj . 

iei iei 

The random variables {/j}j e / are positively related if for each i, there exist 
random variables { Jj^j^i defined on the same probability space such that 

C(J j ,i,j€l)=C(I j j€l\I i = l) 
and, for all j ^ i, Jj^ > Ij. 

Proposition 4.4. If the random variables {/j}j g / are positively related, 
then 

d TV (C(Z),V(9)) < (Var (Z) - 9 + 2 . 

Proposition 4.4 can be applied to our context. First, observe there is a 
natural partial ordering on the configuration set X n = {— ,+} Vn defined by 
a < a' if o~(x) < cr'(x) for all vertices x 6 V n . A function /: X n — > R is in- 
creasing if f(a) < f(a') whenever a < a' . By construction, the indicators 7^ 
are increasing functions. Furthermore, for a positive value of the pair poten- 
tial b(n), the Gibbs measure /^j, defined by (2) satisfies the FKG inequality, 
that is, 

(22) E a , b [fg]>E a)b [f]E ajb \g], 

for all increasing functions / and g on X n \ see, for instance, Section 3 of 
[14]. Then, Theorem 2.G, page 29 of [2] implies that the increasing ran- 
dom indicators I x , x S V n , are positively related. Replacing Ii with T],, 
Z with X n {rf) and 9 with A n , Proposition 4.4 produces the following re- 
sult. 



22 D. COUPIER 

LEMMA 4.5. The following inequality holds: 
d TV {C(X n (r])),V(X n )) 

(23) 

< - 1 - ( Var^^Cfy)] - A n + 2 ^ E a , b [T x } 2 



The last term of the above bound is equal to the ratio A,„ divided by n . 

At, 



E E «, 6 [/ 



■T?i2 



x€V n 

So, it tends to as n tends to infinity. Hence, the bound (23) indicates that, 
as n — > +oo, the distance to the Poisson approximation is essentially the 
difference between the variance and the expectation of X n (rj). This difference 
can be written as M 2 (X n (n)) + A^ where M 2 (X n (rj)) denotes the second 
moment of the random variable X n (rj): 

M 2 (X n (r ] ))=E a:b [X n (7 1 )(X n (r ] ) - 1)]. 

The quantity A^ is controlled by Lemma 4.2, and the second moment X n (r]) 
by the following result. 

Lemma 4.6. The second moment of the random variable X n {rf) satisfies 

M 2 (X n (v)) = A 2 + 0(M n { v )). 
We deduce from Lemma 4.5 that 

d TV (C(X n (r,)),V(X n )) < y (M 2 (X n (r,)) + X 2 J^ - 1 

Thus, combining inequalities given by Lemmas 4.2 and 4.6, it follows that 

d TV (C(X n ( V )),V(X n )) = 0(M n (rj) +n~ d ). 

We conclude using \n~ d = W n {rf) < A n (rj) < M n (n). 

Let us finish the proof of Theorem 3.1 by proving Lemma 4.6. 

Proof of Lemma 4.5. The variable X n {rf) counts the number of copies 
in the graph G n of local configurations belonging to C r {rj). So, the quantity 
M 2 (X n (r])) can be interpreted as the expected number of ordered couples 
of copies of elements of C r (rf): 

u\ 

M 2 (X n (r,)) = Ma,6(*nfa) = k )jT~^\ 
k>2 [K-Z). 



EaJ 



E 



m 



POISSON APPROXIMATIONS FOR THE ISING MODEL 



23 



Let (x\,X2) be a couple of vertices. Prom now on, two cases must be dis- 
tinguished: either the distance between x\ and X2 is smaller than 2r + 3 or 
not. The meaning of this distinction will be revealed later. However, we can 
already remark that the number of couples (2:1,22) such that dist(xi,X2) < 
2r + 3 is bounded by /3(2r + 3)n d . Indeed, there are n d possibilities for the 
first vertex x\ and no more than (3{2r + 3) possibilities for the second one 
since it belongs to the ball B(x\,2r + 3). As for the number of couples 
(xi, x-i) such that dist(xi, x?) > 2r + 3, it is merely bounded by n 2d . 

Each indicator 7^ is defined as the sum of Zj? , rf E C r {rf). So, the second 
moment M2(X n (rj)) can be expressed as 



E 



E ^ x m 



+ E-I 



E 

(zi,z 2 )ev; 2 



&Cr{v) dist(xi,x 2 )<2r+3 dist(xi ,x 2 )>2r+3 

where the above expectations will be respectively denoted by .E- (771, 772) and 
E > (771, 772). We are going to prove the three following statements from which 
Lemma 4.6 follows: 

(24) y( m , m )eC r (r,) 2 E^( m , m ) = 0(M n (r,)); 

(25) E>( V , V ) = \ 2 + 0{M n ( V )); 

(26) V(7 /1 ,7 /2 )GC r (7 ? ) 2 \{(7 / ,7 / )} E>( m , m ) = 0(M n (rj)). 

Let (xi, X2) be a couple of vertices such that dist(xi, X2) < 2r + 3 and let 
7/1, 7/2 be two local configurations of C r {rf). The balls B(x\,r) and B(x2, r) are 
both included in the large ball B(x±, 3r+3). The intuition is that if 771 and 772 
occur on B{x\,r) and B(x2,r), then locally [i.e., in the ball B(x±, 3r + 3)] at 
least k(rj) + 1 positive vertices are present. This has a vanishing probability. 
So as to lighten formulas, let us denote by B the ball B(x\,r) and by R the 
set of vertices y such that r < dist(xi,y) < 3r + 3. The event = 1^ = 1 
implies 



f U U 4=ig=iW u 



where a represents the null configuration of {— First, remark that 
upper bound of Lemma 2.4 can be extended from {— , +} SB to {— , +} R . For 
all 77' G C r (rf) and a G {-, +} R , 

(27) fi a , b (4 = l\a)< W n ( V ') (l a=CT o + t a ^o A '> 



W n {a) 
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Let us consider rf S C*{rf) and a G {— , +} fi . Relation (27) implies 

00,6(4 = = 1) = HaM = MlR = l)Ma,6(JS = X ) 

^-^( 1 + t^)^ (JS = 1) 

where C{R) is the constant of (H2) associated to the set R. Thus, using 

A„(V) <exp(2a(ra)) < 1 
and 77' € C*(rj), it follows that 

A*o,6(4 = = 1) < W n (r])(l + C(R))e n ( V ). 
Now, if the configuration a belonging to {— ,+} R is different from a , then 

VaM = I R = 1) = M J B = ^ = l)Mo,6(JS = 1) 

<W n (r ? )^M0 a>6 (JS = l) 
H/ n (crj 

<W n (r/)C(^)A n (r ? ). 

In conclusion, we obtain an explicit bound for E-(r]i, 772): 

EHm,m) < /9(2r + 3)A(|C;( ?? )|2l R l (1 + C(i?))e n (r7) + 2l fi IC( J R)A n (r / )), 

from which relation (24) follows. 

Let (xi, X2) be a couple of vertices such that dist(xi, X2) > 2r + 3 and let 
7/1, 7/2 be two local configurations of C T (rj). Here, denote respectively by B\ 
and B2 the balls B{x\,r) and B(x2,r). Since the distance between xi and 
X2 is larger than 2r + 3, no vertex of B\ can be a neighbor of a vertex of 
B2 [actually, dist(xi,22) > 2r + 1 suffices]. The Gibbs measure fi a ,b yields a 
Markov random field with respect to neighborhoods defined in (1) (see, e.g., 
[17], Lemma 3, page 7). As a consequence, 

Ma,ft(i£ = I? 2 = 1) = E 0)6 [M , 6 (i£ = ^ 2 2 = 11^(^1 U 552))] 



(28) = E 



0,6 



Li=l 



■ 2 

= Ea, 6 ] 
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In a first time, assume that 771 = 772 — V- Let a be a configuration of 
Thanks to Lemma 2.4, we can write 



8=1 



< W n (r/) 2 ( 1 + A n (r?) ( -— + 



W n (a 5Bl ) W n (a S B 2 ) 

+ A " (??) '{w n (a SBl )W n (<TsB 2 : 

The reason for which the distance between vertices x\ and X2 is assumed 
larger than 2r + 3 is the following: no vertex of 5B\ can be a neighbor of 
a vertex of <51?2- So, the product W n (asB 1 )W n (asB 2 ) 1S equal to the weight 
W n (asB 1 u5B 2 )- Hence, the inequality A n (r/) < 1 and relation (28) imply that 
the probability /U 0) b(7J = £J? = 1) is bounded by 



(1 + A n (r?) (-— — — + — - ^— + -— -i r)/v^) 

Let C(SB) and C(25B) be the constants of (H2) respectively associated to 
the sets 5B\ (or 8B2) and <5i?i U 5i?2- Then, a new upper bound for the 
quantity n a ,b{Ixi = ^x 2 = -0 * s obtained: 

W n (i]) 2 {l + A n (r ? )(2l' 5Bl IC(<5 J B) + 2l <5B2 IC(5 J B) + 2^ 8B ^ 5B ^C{25B))). 

Finally, the term E > (77, 77) satisfies the inequality 

£>M) < n M W n (77) 2 (l + 0(A n (r/))) 

and (25) is proved. 

Assume at least one of the two local configurations 771 , 772 € C r (77) is differ- 
ent from 77. It remains to prove that .£^(771,772) is of order M n (r]). Techniques 
used in the previous case allow us to write 
2 

IIM^ = M^m)) = o(w n ( m )w n ( m )). 
1=1 

Since this upper bound does not depend on the configuration on 8B\ U 8B2, 
we deduce from (28): 

= n\ = 1) = o(w n ( m )w n ( m )). 

Now, if the local configuration 77^, for i = 1,2, is different from 77, then its 
weight W n (rji) is smaller than W n (r])@ n (ri). So, Q n (v) < 1 for n large enough 
and (771,772) 7^ (77,77) imply that the product W n (r]i)W n (r]2) is bounded by 
W n (r]) 2 @ n (r]) and consequently (26) is proved. □ 
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